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' Abstract 

In this note we investigate the impact of new physics in the form of general four-fermion 
{^X)! interactions on neutrino osciUation signals. We develop a field theoretic description of the 

^ . overall oscillation process which includes the non-standard interactions during the neutrino 

' production, propagation and detection stages. Insights gained during the development of this 

formalism regarding the possibility of new interactions mimicking neutrino mass differences are 
expounded. The impact of possible new physics is assessed by studying the i/^ Vr oscillation 
channel in vacuum and the — > v^i, channel in matter. Although it is known that the effects 
of new interactions can only act as a perturbation to the leading oscillation parameters, we 
find that great care needs to be taken when drawing conclusions regarding the strength of new 
couplings from oscillation measurements. 

00 ' 1 Introduction 

O 

' The observation of neutrino oscillations have had a profound impact on our understanding of 

<^ . particle physics. The measurement of atmospheric neutrino oscillations at the Superkamiokande 

' (SK) experiment indicates that f^,'s oscillates to i^r's with a mass squared difference of « 3 x 

10 ^eV'^ , a result confirmed by K2K j2|. The simplest interpretation of this result is that 
at least one of the neutrino mass eigenstates has a non-zero mass of m > 5.5 x 10~^eV, and as 
D . such provides the first glimpse of physics beyond the Standard Model (SM). Recent measurements 

' of the solar neutrino spectrum at the Sudbury Neutrino Observatory (SNO) have confirmed the 

solar neutrino results of SK and found that from the sun are oscillating to the — Vr sub- 
system 13 El El- The corresponding mass squared difference is smaller. 

The standard treatment of neutrino mixing has three active flavour eigenstates mixed with 
' three mass eigenstates. This mixing is described by the rotation of one basis into the other by the 

angles ^13, 9i2 and 623 and a complex phase. The atmospheric data constrains sin023 to be almost 
maximal with a mass squared difference of |<5mi3| « 3 x lO^'^eV^, while the solar neutrino data also 
constrains 6*12 to be large and has |<5rni2| < 3 x 10^'^eV^ 6 . These "large mixing angle results" 
for the solar neutrino mixing were recently confirmed by the terrestrial Kamland experiment j^. 
Furthermore 6*13 is bound by reactor data such that sin^ ^13 < 0.1 0. 

The signs of the mass differences are still unknown, while information on possible complex CP 
violating phases is extremely scant. In order to overcome this there have been a number new 
experimental facilities proposed such as the the JHF neutrino beam, and others based on a muon 
storage ring. All provide a high intensity neutrino beam over a wide range of base lengths and 
energies |3[TnilIIl. 

Alongside the neutrino oscillation industry there has been an ongoing and long established 
program of precision experiments aimed at measuring the properties of the weak interaction. Mea- 
surement of the Tritium beta decay spectrum end-point is one such program while the determina- 
tion of the Michel parameters in the /i and r-decay spectra is another |12II13I .14, 15, 16, 17 1^. 
Originally the focus of these experiments was to establish the Lorentz structure of the weak inter- 
action, that is to decide whether it is Vector minus Axial- Vector (V-A) or some other combination 
of Scalar and Tensor operators ^^1- Now that the V-A structure has been established the focus 
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has shifted to measuring the electron neutrino mass, at least in the case of the Tritium experi- 
ments 1201. Not until very lately have the possible effects of neutrino mixing on these experiments 
been considered [HI 122123 Ell- 
in this paper the converse scenario is examined, that is the impact of non-standard interac- 
tions, such as additional Lorentz structures in the weak interaction, on future neutrino oscillation 
experiments. The aim of this work is to derive a formalism that will describe the production, 
propagation and detection of a massive neutrino with general Lorentz structures present at each 
of these stages. In effect we aim to incorporate the couplings that comprise the Michel parame- 
ters, and the /3-decay parameters of Ref. 19 into the description of neutrino oscillations. Similar 
calculations have focussed on specific extensions to the Standard Model, or have been performed 
exclusive of one or more of the production, propagation or detection process. The correspondence 
between these previous calculations and this work is expounded. 

2 Background to NSI 

In the SM weak processes are mediated via the exchange of a charged vector boson, and are 
successfully described at low energies hy a V — A current-current interaction with an effective 
coupling strength Gp- Presently the SM has withstood all experimental tests with the exception 
of the recently observed apparent non-zero mass of the neutrino. Despite this there are many 
theories that seek to extend the SM, often motivated by the desire to restore a broken symmetry 
of the SM or to unify the strong force with the electroweak force, or even to unite all four of the 
fundamental forces under one theory. Usually this unification occurs at a higher mass scale than 
the SM, that is the bosons that mediate the new interactions have a mass larger than the W and Z 
bosons of the electroweak theory. The new interactions will manifest themselves as current-current 
theories at low energies in the same way as the SM charged current reduces to Fermi's effective 
field theory for energies significantly less than the mass of the Z-boson. 

If new physics is present at higher mass scales the low energy result will be to induce new 
effective interactions with a coupling strength weaker than the coupling of the dominant V — A 
interaction unless the fundamental coupling constant is anomolously large. Absent a specific model 
there is no reason to expect that the new physics has the same Lorentz structure as the weak 
interaction, hence new currents may be a scalar (S), pseudo-scalar (P) or tensor (T) in nature or 
comprise of different combinations of vector and axial vector operators. Given the dominance of 
the left-handed V — A interaction we find it convenient to cast these new interactions into forms 
with definite handedness\ L = V - A, R ^ V + A, Sl = S - P and Sb. = S + P. 

Assuming lepton universality the most general interaction Hamiltonians for low energy, leptonic 
and semi-leptonic processes are given by 

Hl = ^g^'^V^iP/sV^^^ (v;™r„V,o.0)V/i.c., (1) 

a/3 

and, 

HsL = ^G^'^^iP^V,"-^ (VigPa^^O^ + Z^.C. (2) 
a/3 

respectively. Where a//3 — L, R, Sl, Sr, while ipq and ipgi are quark fields, -ipi and ipm are 
charged lepton fields of flavour / and m respectively. The neutrino field is a^, representing a 

neutrino produced by an {a(3) type interaction and associated with a charged lepton of flavour I. 
The operators Pa are a combination of one of the flve bilinear covariants 

Pl = 7^(1 -75) , 

Pfl = 7^(1+75) , 

Ts. = (I-75), 

rs„ = (I + 75). (3) 

^To be completely general one would need to also consider new V-A interactions in addition to the standard 
interaction. This has been the subject of numerous investigations, see Ref 1251 1261 1^71 . 
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In the SM only a — {L,R) and [3 — L are present in the semi-leptonic case with tl^u^^ = V'l/"^- 
And for the standard Icptonic interaction only a — L and (3 = L are allowed. 

The fact that neutrinos oscillate indicates that the interaction eigenstatc is not the same as the 
mass eigenstate. Typically the two are related via a unitary transformation 

i^.,^Y.Uf.i>u. ■ (4) 

i 

Phenomenologically, there is no reason why a neutrino produced by one interaction needs to be 
coupled to the same combination of mass eigenstates as a neutrino produced by another interaction. 
So in this work neutrino mixing is described by 

i 

Or equivalently each new interaction has its own unitary mixing matrix which may be equivalent, 
or not, to the standard mixing matrix. 

The Eqs.n& El describe the production and detection of a neutrino associated with a charged 
lepton, in addition Eq. ^can describe the interaction of a neutrino, via a charged current, with 
a background medium of charged leptons. This interaction leads to the well known resonant 
enhancement of the neutrino mass, the MSW effect, crucial to the solution of the solar neutrino 
problem (2S1 129| . In order to be used in the derivation of the matter induced potential, to which 
we will return in a later section, Eq. ^ needs to be cast in the form of a neutral current. This is 
achieved by Fierz rearrangement (See for example Ref. (2111) such that 



1/3 
a,0 



where 



a,f3 ij 



n-RR _ ^RRttRRttRR 



r^RL ^ ^SlSlttSlSlttSlSl 

- 2^ '^u '^ij ' 



^m^ = ^.9^''-^^ ut^"- Ui,''^'- + 35^^-^^ ul^^^ Uj'^^^ , 

T^SrSl r\ RLttRLttRL 

~ ^ ' 
t^SlSr c\ LRjjLRjjLR 

t^TlTl ^ „Sr,Sl ttSrSl ttSrSl ^ „TrTl ttTrTl ttTrTl 

(Oy ~ 4^ iJ 2 '-J' 

T^TrTr ^ „SlSrttSlSrttSlSr ^ TlTrttTlTrttTlTr /i^N 

(Oy ~ y 2 '-J' ■ ^ ' 

This potential has been written in the mass basis so that Eq.Elactually describes the interaction 
of a mass eigenstate, rather than a flavour eigenstate, as is usually presented. The effective coupling 
constant K'^^^^. is dependent on the type of leptons in the background as it contains the various 
mixing elements specific to the lepton fiavor. We do not included a neutral current term for two 
reasons. Firstly, the standard neutral current contributes only to the absolute value of the effective 
neutrino masses and not to the mass differences or mixing angles. And secondly, any non-standard 
neutral current is highly constrained by experimental data from atomic physics and LEP. 
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In this work we apply the formahsm being developed presently to terrestrial experiments, 
hence a background of charged fermions only is considered. In the general case one would also 
need to consider the effects of forward scattering off a neutrino background. This has important 
consequences in early universe physics, described for example in ''M\ 1^211!^ . 

As alluded to earlier the Lorentz structure of the various currents in the effective Lagrangians 
presented above may be generated by the presence of physics beyond the SM. One prominent 
example is the low energy effects of Supersymmetry |34| . In particular R-parity violating modes 
lead to charge-changing scalar currents 35 . In the model discussed in Ref. I^HI stau mediated 
lepton number violating decay, /i^ — > e~i^eP^, is examined. For instance the correspondence 
between the Lagrangian that produces this decay and Eq. ^may be found by setting ^^s^s^ = 
tp^LL and '0 ^L^L = V'l/J^i and if for the supersymmetric coupling the following substitution is 

made Ai32A23i/(s — m^) = G^^^'-, where is the stau mass and A132 and A231 are coupling 
constants. Supersymmetry aside, right-handed vector currents can be motivated by left-right 
symmetric models where a heavier W boson couples to right-handed neutrinos, while scalar, vector 
and tensor currents may arise from Leptoquark theories j36| . 



3 Formalism — Scattering theory 

Spectral measurements such as the Tritium and Michel parameter experiments yield little infor- 
mation about the mixing elements U'^f , hence it is natural to look to oscillation experiments to 
access these parameters. To do this a formalism that describes the oscillation of a neutrino with 
information about both the production and detection processes is needed. To do this we treat the 
whole production, propagation and detection process as a single scattering event. In this scenario 
the neutrino and charged vector bosons at the production and detection sites are treated as un- 
observed intermediate states. The oscillation phenomena arises through the interference between 
these scattering diagrams. Now if NSI are present additional diagrams with the non-standard 
boson at the production and/or detection site need to be included. 

Field theoretic (scattering theory) calculations of neutrino oscillations have been performed by 
a number of authors [371 OHl EHI EHl E] ■ For a complete list of references and a thorough review see 
the work by Beuthe, Ref. In this note the formalism developed by Cardall & Chung, Ref. [31] 
is adapted to cater for NSI. 

To keep the notation transparent the formalism will be developed for a specific example, that 
of /i^-decay at the source, electron neutrino^ oscillation and subsequent detection of a negatively 
charge muon at the detector. The scattering amplitude is given by the time ordered product of 
interaction Hamiltonians 



A = <i?^e+Xti-\T[J d^x J d''yHLix)HsLiy)]\tJ-'^N> 
dx° / dy" / d^x / d^y 

- 00 J —OQ J V.q J Vr> 



^ Vi2E,.)i2E.M2E,) ^/'(--'^l^-^/-'^) ^i2Ex)i2E,-)i2Ex) ' 

where Ee and Ep, are the energies of the anti-muon, electron and the muon anti-neutrino 

(not participating in the oscillation measurement) with corresponding four momenta , Pe and 
Pp. The reaction at the source is evaluated at a space-time point x in a volume Vs- The energies 
of the particles at the detector are E^, E^- , Ex with associated four-momenta pat, p^^- , px- The 
reaction at the detector is evaluated at a space-time point y in a volume Vd- Furthermore the 
matrix element Mfi{x,y\pi,pf,q) is given by 

Mj,ix,y\p.,Pf,q) = J2J2J29"''G^''KfU^IJ^iP^^P^')^ 

ij a.(3 Act 

X {u{p^ -)TM'\x,y)fTifv{p,)) ■ {v{p^+)T0v{p,)) . (9) 

^The term electron neutrino is used loosely here since neutrinos produced via different interactions and associated 
with the decay of muons are are not necessarily the same. This is due to the various mixing matrices associated 
with the different interactions. When it matters more specific language will be used. 
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Here pi and pf are the momenta of the particles in the initial and final state respectively, q is 
the momentum of the oscillating neutrino while J^{pN,px) is the nuclear matrix element for 
the transition operator T\. The sums are over neutrino vacuum masses, source and detector 
interactions respectively. The function G*-' (x, y) is the neutrino propagator or Green's function. In 
vacuum, G^^{x,y) G^'^{x,y) and is given by the standard Dirac propagator. 

The starting point in the evaluation of G'^ (x, y) is the equation of motion for a propagating 
neutrino from which the equation satisfied by the neutrino Green's function can be found, 

{i^-M-V)Giy,x)^S^iy-x) , (10) 

where V is the propagation potential derived from Eq. and M is the mass matrix, note the mass 
indices have been suppressed as we will be concerned with G(y, x) in spinor space for the moment. 
Due to the chiral nature of the operators and Fj^ it is convenient to write the Green's function 
in chiral blocks 

G(x,y)=f ^^^[^'^J ^^^["'^l ) , (11) 
^ \ GRL{x,y) GRR{x,y) J ' ^ ^ 

where Gxy is the element projected out by PxG{x, y)PY where Px — Pr/l ~ Similarly 
the operator in Eq. IIUI in chiral form can be written as 

where 

Ml = M + S*^"^^ (13) 



and 



with 



Ma = Af + , (14) 



The sums in Eas. 115111^ are over lepton species with mass mi distributed according to pi{pi) in the 
background medium. In general the Green's function may be expressed in any basis, however we 
find it convenient to do so in the basis defined by the vacuum mass eigenstates. In this basis the 
matrix M is diagonal, while the matrices S^"-^^, S^'-^"-, (y^^)^ and (y^^)^ are in general not 
diagonal. All other operators are proportional to the unit in this basis assuming the potential V 
does not vary in space or time. 

Using Eq. & Eq. El two sets of two equations can be found and used to solve for each 
component of the Green's fmiction. For Glr and Grr; 

S\x2-xi) = -AfsGi^fi + [z(9o - • V) - • 5] Gla , (19) 
= -MiGLR+\iido + cj-V)-V'^''-a\GRR, (20) 
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and for Gll and Grl] 



5^{X2 - Xi) 







i(5o - 
i(9o 



CT- V) 
- • V) 



Grl 

Gll 



(21) 
(22) 



For the analysis that is to follow in the next section only the first set of equations will need to be 
solved. This is done following Cardall & Chung, Ref. by defining 



Jix2,Xi) = Afj ^Grr{x2,Xi), 

which reduces the problem to solving 

(-M2M1 + [ia.d -VL-a][ia-d-VR-a]) J{x, y) = (5*(y - x) 



(23) 



(24) 



where vr, = M-^ ^V^^Mi and for brevity vl — V^^. Now the delta function and J{y,x) are 
expanded into momentum space via a Fourier transform where we find that 



{-M2M1 + [{qo + (T ■ q) - VL ■ cr] [(go - S ■ q) - vr - (T])j{qn,q) = 1 



(25) 



Where qo and q are the unobserved neutrino energy and momenta. The task is now to find the 
elements of J{qo,q} in terms of its reciprocal matrix R{qo,q) = J{qo,q}~^. Recalling that J{qo,q) 
is a matrix in the 2x2 space of the Pauli spinors, one can make a substantial simplification by 
choosing the direction of neutrino propagation to coincide with third spatial coordinate. For a 
non-polarised, non-relativistic background of leptons we find that: 



\q\)i2 
R{qo, \q\)2i 
R{qo, \q\)22 



ql 



qovl-vl\q\]+vlvl, 



ql - |gf - M2M1 ~ [go< - 

[qovl + vim+vlvl, 



(26) 
(27) 
(28) 

(29) 



where the subscripts correspond to the elements of the matrix R{qo, q)- 

To perform the Fourier transform back to coordinate space the pole structure of J{qQ, q) must be 
established. This is accomplished by defining the unitary matrices Ur and U l which diagonalize, 
in flavour space, Ru and R22 respectively. In the ultra-relativistic limit where \q\ qo, the 
components of J{x2,xi) are 



J{y,x)22 



dqo 



-tqaiy^-x") 



"■RK I.-? j;l 



\v- 



2(27r)2 
dqo 



\y - x\ 



2(27r)2 



\y - x\ 



K 



(30) 
(31) 



K 



The effective masses mj^j^ and mj^j^ are the K eigenvalues of the matrices —M2Mi+v'^m'~' — 2qom^ 
and — M2M1 + w^m*^ — 2(jow° respectively. 

The components of the Green's function, Grr and Glr, are found by substitution into Eq.^S 
and then Eg. 1191 In the matrix element, Eq. [SJ the sums over vacuum masses are carried out such 
that 



(32) 



K 



where the sum is now over effective masses. 
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The scattering amplitude is found by localising the particles at the source and detector to boxes 
of volume Vs and Vd and by expanding x and y about the source and detection points xs and yjj 
respectively PII- Finally we find for the square of the scattering amplitude, 

Vs Vd 

X {27r)'^5^{p^+ -Pe~Py- q){2TTf6'^{pN + q - p^- - px) 
/ \ 2 

exp[z5m,2^|L|](A^^+A^f;)(X^_Al^_*) . (33) 



To find a per source particle per detector particle event rate integrate over the final state phase 
space (factor of Vd?'p/{2'KY' for each particle in the final state) and divide by T after interpreting 
one of the energy delta functions as characterising the time of the scattering process; 



|i|2 dEqdVLq 



KL l-l 



a{uN)KLdEq , (34) 

KL 



The quantity d'^N^/dEqdO.q is interpreted as the neutrino flux; 



jqu^uq 

d^N,. 



dEqdQ.q 



KL 



d^p, d^p^ {'27r)HHp^+ ~Pe-p,- q)E^qM^^tM^ 
(2^)3 (27r)3 {27Tf{2E^^){2Eq){2Ei.)(2E,) 



K 
+ 



(35) 



where the {KL) dependence arises from the various mixing elements in matter. The {vN) cross 
section also has a [KL) dependence and is given by 

I m f '^'P^ '^'P^- i^^)'^HPN + q-Px-P,-)M^^-M;^* 

~ JWrlW i2Eq)i2E,)i2Ex){2E,-) ' ^''^ 

The mixing elements are all contained in the squared matrix elements and are not in general able 
to be extracted and interpreted as an oscillation probability, except in the case where only one 
type of interaction is responsible for source and detection processes, or if only one type of mixing 
matrix exists, i.e. U"^ = . 

The event rate is obtained by integrating over the source and detector particle distributions: 

dE-q^J'^'^'J T^/-^^'^-'^^) J^y^j j^f^'^p^'y^Wq ■ 

This event rate, with Eas. 1201 and 1221 is the first main result of this paper. It is an expression 
incorporating neutrino oscillations in matter, and after following the same procedure just outlined 
to solve for Glr & Gll, allows for arbitrary Lorentz couplings. There are three main points of 
difference between this result and the standard result. Firstly, in vacuum there is the possibility 
of observing flavour violation as a result of the intermediate neutrino even if the neutrino mass 
differences are zero. This requires that the non-standard mixing matrix not be related trivially 
to the standard matrix. In this scenario no oscillatory phase will develop. Secondly, in matter, 
the addition of NSI allows for an oscillatory phase even if the neutrino masses vanish. Again the 
non-standard mixing matrix must not be trivially related to the standard matrix. These results 
were first noted by Bergman et. al. in Ref. 42 and subsequently investigated by Huber et. al. in 
Ref. 03! • They found that NSI could only act as a perturbation to the solution of the solar and 
atmospheric problems rather than a complete explanation. And thirdly, we note that the presence 
of a right-handed potential connects a left-handed neutrino state to a right handed state, even 
in the ultra-relativistic limit. This is akin to the spin flip caused by a non-zero neutrino magnetic 
moment in a strong field. We do not pursue this here, however we note that the presence of this 
potential in a dense medium will require the equations of motion for right-handed and left-handed 
neutrinos be solved simultaneously, perhaps using a density matrix approach. Finally we highlight 
the fact that if the right-handed potential vu is zero then the only non-zero component of the 
Greens function is G22 in the ultra-relativistic limit. Thus the connection to right-handed states 
vanishes and the only required diagonalization matrix is . 
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4 Considerations for a neutrino beam 



The formalism developed in the previous section can now be applied to a generic neutrino factory 
scenario. In this section no attempt is made at a full simulation of this type of experiment, A 
full simulation would involve a discussion of detector properties, efficiencies, energy resolutions, 
cuts, backgrounds and the like; rather we attempt to define regions of interesting parameter space 
and achieve a quantitative understanding of the effects of NSI. As we have seen, in the relativis- 
tic limit, the event rate for a particular experiment can be obtained by evaluating the neutrino 
flux, (PN^/dEqdil.g\KL, containing the effects of NSI at the source, the production cross section, 
a{vn)KL, accounting for NSI at the detector, and the phase factor exp(i(5m|-^|Z|/2_Eg)/|ip deal- 
ing with the oscillatory phase. In this paper we present the results of two investigations. Firstly 
the vacuum propagation of a neutrino produced at the source via /i~ -decay and subsequent pro- 
duction of a T~ at the detector, that is z^^ — *■ oscillation. Secondly we examine the case of a 
neutrino produced via /^^-decay, interacting with the electrons in the background medium, and 
producing a /i~ in the detector, or v,, — s- v^^ oscillation. In both cases this is performed for the 
simple case of one non-standard scalar coupling, {a, (3) = {Sl,Sl), in addition to the standard 
interaction in Eqs.^&IJl This particular non-standard coupling is chosen since it is the only inter- 
action other than the standard one that does not require the existence of right-handed neutrinos. 
The strength of the coupling for both the leptonic and semi-leptonic interaction is taken to be 
gSLSL — G'^iS'i. — 0.01 G^^, within the current experimental upper bounds 44, 45^. 

The expression for the neutrino flux resulting from the decay of a /i~ can be calculated using 
standard trace techniques, in general one will find two new terms as a result of the scalar interaction, 
an interference term between the standard and non-standard interaction of order |(?^^'^^| and a 
pure scalar term of order \g^^^^\'^. The interference term refiects the right-handed admixture, 
proportional to in the left-handed wavefunction of the ^~ 0^1 • The chiral structure of the 

scalar interaction being considered requires a right-handed /i~ , while the standard weak interaction 
always proceeds via a left-handed pL~ . Since the process by which muons will be produced at a 
neutrino factory will be dominated by the standard weak interaction, e.g. 7r-decay, any contribution 
from NSI will be doubley suppressed, once by the production process and once by the decay process. 
For this reason we neglect the contribution from NSI at the source. However the expression for the 
neutrino flux with the additional NSI terms is recorded here for completeness. In the rest frame 
of the muon, neglecting all lepton masses we find that 



where x = ^m^j^jE^ and p = G^^/G'^^'^^ . 

For energies above the tau threshold of ~ 5GeV the dominant reaction mechanism at the 
detector is deep inelastic scattering producing many particles in the final state in addition to the 
tau, denoted collectively as X. The {vrN) cross section can be derived using a parton model of 
the nucleon. In the standard treatment the lepton and parton masses are typically assumed to be 
zero, this is a good approximation in most circumstances however its validity for this application 
is questionable. This is due in part to the considerable mass of the r and also to the fact that 
certain mixing patterns may act to pick out the interference term in the cross section which would 
otherwise vanish in the massless limit as we shall show. A treatment using massive partons has 
been given by Aivazis et. al. in Ref '47' , and the non-standard cross section is derived using this 
formalism. 

The general expression for the cross section with a vanishing neutrino mass is 



where A(s,0,M^) is a kinematic function of the Mandelstam variable s and the nucleon mass 
M, while we have retained a non-zero T-mass we have neglected the small neutrino mass. The 



d^N, 
dxdflq 



(38) 



da 




(39) 
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final state phase space is represented by dT. The complete evaluation of Eq. (|39|l is lengthy but 
involves only standard trace algebra, here we just state the results. The cross section in terms of 
the Bjorken scaling variables x and y is written as 



da 



dxdy 



jjLLjjLL 



da 



dxdy 



LL 



ttSlSl ttSlSl 
rj ri 



da 



dxdy 



ss 



_ I tjLLtjSlSl 



tjSlSl TT^L 
rj ^ ri 



da 



dxdy 



(40) 



SL 



where the first two terms arise from vector and scalar interactions respectively and the final term 
arises from the interference between the two. In the limit where all terms of order mf /Q^ or greater 
can be ignored (where in this context rrii corresponds to the parton masses) we find 



da 
dxdy 



\G 



SlSl 



SS 



27r 



xy 



mly 
2ME^ 



Fi{x), 



(41) 



da 



dxdy 



SL 



qS,s,qLlME^ 

277 



xy 



E,. 2ME?/ 



Fiix), 



(42) 



The term, [da/dxdy]^^, not shown here is the standard results and was derived in Ref 0S|- The 
structure functions, Fi through F^, are measured quantities fitted to a functional form, we use the 
same form as is used in Ref ^H]- In the limit that ra\/Q^ the struck parton mass becomes 
mi xM as noted in Ref j5()| . 

With the assumption of an idealised source of neutrinos arising from muon decay the charged 
current event rate can be written as 



Nr 



6.023 X 10^2 



N.Mkt 



E 



d^Pt 



dxdVlr, 



ayN\ij exp(-i- 



2Er, 



dE, 



(43) 



where uj^d^p*/dxdflq — d^N* /dxdflq is the neutrino distribution in the reference frame of the lab. 
For high energy muons x Eq/E^ and: 



dxdflf, 



d'^Pu 



7-^(1 — /3 cosa)2 dxdQ.q 



47- 



dxdQ 



(44) 



where a = coincides with the beam direction, 7 = E^/m^ and (3 — p^/E^. For high beam 
energies cos a ^ 1 and [3^1 — 1/(27). The number of useful muon decays per year, iV^, is defined 
as the integration over the source distribution times the /i-decay rate. The parameter M^t is the 
mass of the detector in kilotons. The numerical factor in Eq. (|43(l is the number of nucleons in the 
detector per kiloton, we have assumed an ideal detector comprising 6.023 x lO'^^Mfet non-relativistic 
nucleons. 



5 Vacuum oscillations — Two neutrinos. 

We now study the effects of varying the non-standard mixing angle for the — > v-r vacuum 
oscillation channel. 

The study is based on the philosophy that the atmospheric, solar and Kamland neutrino data 
sets have accurately defined the leading oscillation parameters. That is the LMA solution is 
accurate and the NSI are treated as perturbations to this solution. This philosophy is backed by 
the study of atmospheric neutrinos and NSI in Ref [SI]. We examine a simplified two neutrino 



9 



system through the i/^ —f v-t osciUation channel. As such the standard and non-standard mixing 
matrices are defined as 

C/^^f '"^M and U^-^-^i ""^'^^ '^^t'l . (45) 

\^ — sm {) cos B J y — sm q> cos (p J ^ 

The standard mixing angle and mass difference given by the LMA solution are sin^(26') « 1 and 
« 2.5 X 10-3 eV^ 

The baselines over which it is proposed that oscillation experiments will be conducted vary from 
a few hundred kilometres to several thousand. For example at the proposed Fermilab neutrino 
factory there are plans to place detectors 732 km away at the Soudan detector and even at the 
South Pole a distance of 11700 km. The energies of the muon beam used to produce the neutrino 
flux are to be optimised for the detection of CP violation in the neutrino sector. Typically energies 
of 20 — 50 GeV are being studied. The great advantage using a muon storage ring as a neutrino 
source is the extremely high intensity of the resulting beam with ^ 10^^ neutrinos expected to be 
produced per year. 

Other studies have investigated the effects on oscillations of varying the base length |L| and the 
beam energy Eq with NSI present |25j . In any case variations in the beam length will have little 
consequence for the parameter range we will consider; that is for energies, mass difference and base 
lengths where 5m|^- |Z|/2£'g << 1. For the v^^ — > oscillation channel in the absence of NSI 



1 / 



2N 



P^,r ~ 2(cos6'sin6')'^ - 2(cose'sin6')M 1 - 

= (cos0sin0)M^|Uj . (46) 

When this is substituted into an event rate such as Eq. the dependence on \L\ is removed. 

The sensitivity of a future neutrino factory to NSI is investigated by defining a function 
which determines the required detector mass and number of useful muon decays in order to claim 
new physics. The function assuming Gaussian statistics for a detector mass and muon number 
of TVp • Mdet ~ lO^^kt/yr is defined as 

|2 



XNM21 - 2^ 77sM ' K^') 

k k 



where is the expected number of t producing charged current events in energy bin k in the 

absence of NSI and for this detector mass and muon number 

N^" = 6.023 X 10^2 ^^^^^ 



dayN\ij exp(- 



2Eq 



dEg . (48) 



Here the energy bins are defined such that Ek-i < Eg < Ek {k — 1, 2, 3, .., n). Furthermore N^^^ 
is the number of charged current events in the fc*'' energy bin expected with the NSI present. 
We define a required detector mass-muon number unit, NM^-ec = 1 x 10^^ kt/yr, to quantify the 
sensitivity to NSI for a given set of parameters. The constraint on NMrec is 

^M.ec > (49) 

where Xgg% is the value at 90% confidence level and one detector mass-muon number unit defined 
as a function of the number of degrees of freedom'^. The sensitivity of NMr. 

■^The subtleties of performing this kind of analysis have been examined in Ref. 1251 . 



10 



energy bins used, where the number of bins corresponds to the number of degrees of freedom, was 
investigated by varying the number of bins for various values of beam energy. We found the result 
to be statistically stable for ten bins or greater. 

In Fig. n the variation of NMj-ec with respect to the non-standard mixing angle (f> is shown. 
The first plot is for a beam energy of = 50 GeV , while the second is for a beam energy of 

= 20 GeV. Points of particular interest are 4> = 7r/4, 57r/4 and (j) — 37r/4, 7tt/A. The first case 
corresponds to the situation i^^^^^ = or C/^/"^^ = C/^^^ up to a phase. For this case NM^ec is 
close to a maximum. In the second case v^^^^ = v^^ or U^^^'^ = U:^/" up to a phase, these points 
correspond to direct flavour violation. In this case NM^ec is almost at a minimum. This becomes 
obvious when we consider the mixing of the non-standard basis within the standard by performing 
the rotation {v^^) = u^^^'U^^^^ (v^^^^) = I^(z/|^'^^); now 

( ^i"^ \ ^ ( '^'''^'^ ^ '^"'^'^ '^^'^'^ ~ '^"'^'^ ^ ( ^M^'^^ 

\ U^^ ) \ SgC^ - SgS^ + CgC^ J \ v^^Sl 

For the case of sin(26') = 1 the standard basis coincides with the non-standard when cos((/)) — 
sin((/)) = up to a phase and is directly flavour violating when cos(0) -I- sin((/)) — 0. 

The two plots also exhibit a strong energy dependence with a much richer structure evident 
for the higher beam energy. This effect may be understood through a careful examination of the 
{vrN) cross section. The interference term is suppressed by a factor of E~^ relative to the other 
terms, thus its importance is enhanced at low energies. In Fig. |21the relative contributions to the 
event rate from the pure scalar and the interference terms are shown. The flrst plot is with a beam 
energy of E^ = 50 GeV and the second with E^ = 20 GeV. At high energies the magnitude of 
the contribution from the pure scalar term is approximately equivalent to that of the interference 
term, this is despite the fact that it is suppressed by a relative factor of G^'^^'^ . At the lower beam 
energy the interference term is not yet washed out and dominates over the pure scalar term. 

This simplified analysis indicates that for some values of (/) a detector mass of ^ lOOOkt would 
be in a good position to either detect NSI or increase the upper bound on the non-standard 
coupling strength. This result is to be expected, since at this energy and over a medium base 
line the neutrino beam will comprise mostly of non-oscillated t'^'s. Any flavour violating non- 
standard coupling will in essence be picked out over the standard term. This example effectively 
demonstrates the convenience of this formalism, it allows one to obtain results for flavour violating 
couplings, as in Ref. or flavour diagonal couplings simply by dialling up the appropriate value 
of (j). It also serves as a useful reminder that absent a specific SM extension the phase of a new 
interaction can play as important a role as the coupling strength in oscillation experiments, while 
it is unimportant in spectral tests of NSI. 




6 Matter enhanced oscillations — two neutrinos 

We now move on to examine the impact of NSI on the experimental signature of the matter 
enhanced i/g — > i/^ oscillation channel. Again this study is conducted with a generic neutrino 
factory in mind. For this calculation only new interactions during the propagation stage of the 
oscillation process are considered, no new interaction are present at the detector or the source. In 
particular a left-chiral scalar interaction coupled with a strength of g^'-^'- — O.Olg^^ is examined. 
Since there is no new physics at the detector or the source it is possible to factor out an 'oscillation 
probability' in Eq. 021 We will comment further on this issue later. 

The propagation potential for a neutrino with the SM and left-chiral scalar interaction present 
is given by Eq. (TJ and Eq. ()17|l in the mass basis as 

.7^(1-75) = A/2Gj.t/f,^C/i;-V-^^Gj.f/i;-^-C/i;-^-ne (51) 
or in the SM interaction basis 

S S 

C^")^ -7^(1 -75) = V2GFne-y^^-^GFW^,W0,n, (52) 

t^ff V2 
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where W = ijLL-\jjSlSl -^^rg^g introduced in Eq. [SDI Wc have denoted the non-standard basis 
with the tilde in Eq. H52I) . The propagation potential is rotated into the SM basis for ease of 
comparison with the standard treatment of matter enhanced oscillations, in this spirit the matrix 
W is parameterised in an analogous way to the non-standard vacuum mixing matrix: 



W 



cos ip sin ^ 
— sin Ip cos "0 



With these considerations the Hamiltonian for neutrinos in matter is 



H 



2E„ 



U 



LL 



^ 
2 " 





A 




(53) 



(54) 



where A' = 2V2^-4rL^GFneEg and A = 2^/2GFneEq. The expression H can be rewritten as 



H 



^E„ 



A 



A 



A 
~2 



A_ 

2 



5m 



- cos(2^/') snv{2ip) 
sin(2?/') cos(2-!/)) 



~-cos(26') sin(26l) 
sin(26') cos(26') 



(55) 



where E — m\ -\- m\ and as usual 5m? — m% 
basis such that Hm ■ 



if. Now an angle x that rotates H into a diagonal 
Um{x)HU{x)m is found to be 



tan(2x) 



A' sin(2-0) + 25m'^ sin(26') 



25m^ cos(26') -2A + 2A' cos(2?/') 
There are two distinct eigenvalues of Hm 



(56) 



rriLK = ±1 



^ A' 
A — — cos(2V') 



5m^ cos(26i) 



Y sin(2i/^) 



(5m^ sin(26') 



A- 



A 



(57) 



The Eas. l56l and l57l elucidate why oscillations are possible even for vanishing or degenerate vacuum 
masses. A non-trivial phase develops and the mass degeneracy is lifted due to the non-standard 
interaction in the medium, ie the effect is proportional to A . 

The validity of performing an analysis of matter enhanced oscillations with NSI using a two 
neutrino model may be questioned. The approximations that allow the two neutrino analysis to 
be performed in the standard case need not hold in the presence of non-standard interactions. In 
the standard oscillation scenario the three neutrino model can be described with two neutrinos for 
some values of mixings and mass differences. The standard 3x3 mixing matrix is written as the 
product of three rotations U = R{923) ■ ^(6*13) • R{di2)- Typically the mass difference 5m2i is taken 
to vanish approximately. This decouples the R{0i2) rotation, in addition the i?(023) operates in 
the 23-subspace and commutes with the matter induced term, meaning that only ^13 is enhanced 



by the induced potential. For the 



channel the standard oscillation probability becomes: 



Pe;. = sin^ 023sin^(20f|)sin^( 



5m 



13 



iE„ 



(58) 



where df^ is the effective mixing angle in matter and (Jmfg is also the effective mass difference 
in matter. Use can be made of Eq. EHl in studying the effects of NSI with the mass difference 
and effective mixing angle given by Eq. 153 and Eq. 1^01 if one assumes that the NSI basis is only 
non-diagonal in the 13 and 12-subsystems of the standard basis. Furthermore one has to assume 
that the non-standard equivalent of R{0i2) is approximately equal to the identity. For the moment 
we acknowledge these assumptions as potential pitfalls and perform an analysis for the sake of 
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gaining an intuitive understanding of the system. The three neutrino scenario will be examined 
in the next section. For this calculation a standard parameterisation of the {v^j^N) deep inelastic 
scattering cross section may be used since we are assuming no new physics at the detector, that 

is [ig 

9 

cm 

a,^{E,)c^OmxlQ-^^E,-^ . (59) 
Furthermore the electron neutrino energy distribution has the standard form 

=lia.2(l-a;). (60) 



dxdVlq An 

The matter induced potential can be written as 

A = 2V2GFYepEq = 1.52 x lO-'^eY^Yep{g/cm^)Eg{GeV) , (61) 

where Ye is the electron fraction and p is the density of matter along the neutrino trajectory. 
Typically ~ 0.5 in the earth. The value of p varies between ~ 3 — 4.5g/cm'^ depending on the 
base length and tilt angle of the experiment. 

The dependence on ip is examined in Fig. This plot shows two curves one with i?^ = 50 GeV 
and the other with E"^ — 20 GeV. The minimum value of NMrec occurs for tp = tt/A, as one would 
expect since this angle corresponds to a maximal mixing of non-standard flavours in the standard 
basis. In fact the behaviour of the diagonalization angle x for high energy neutrinos also leads one 
to conclude that the greatest effect will be for ip = 7r/4. For large neutrino energy and tp = 7r/4 
Eg. ICT becomes: 

tan(2x) 



-2A 

= -0.005 , (62) 

for a non-standard coupling strength of g^^^^ = O.Olg'^'^. Without the NSI in the higher energy 
limit tan(2x) — > 0. These results are in contrast to the vacuum oscillations studied in the previous 
section. The case of new physics at the detector showed the greatest sensitivity when ipi — {^^^). 

Given that the required detector mass-muon number varies so dramatically with mixing angle 
one must conclude that it is not wise to make the approximations that allowed the study of 
this system in terms of two neutrinos only. In addition, note that the effects of NSI during the 
propagation stage of the oscillation process are about three orders of magnitude greater than for a 
non-standard interaction of the same strength at either the detection or production stages. This 
implies that one can study NSI during propagation in isolation of the detector and the source for the 
matter enhanced channel, an assumption made without justification in previous studies |51[I43[I^ . 



7 Matter enhanced oscillations — three neutrinos 

The cautionary comment of the previous section stating that the assumptions required to treat 
the NSI in a two neutrino framework is investigated by extending the formalism to allow for three 
neutrinos. This is done presently, in addition to quantifying the assertion that for matter enhanced 
oscillations NSI at the source and detector may be neglected. 

Extending the formalism of Sectional presents no real hurdles, but an analytic understanding of 
the results is now difficult. We find it convenient to put the standard and non-standard interactions 
with the lepton background on the same footing by evaluating the propagation potential. Eg. 1511 in 
the mass basis. The effective mass of a neutrino in matter is found by evaluating the eigenvalues of 
the equation of motion in this basis. The unitary matrix, Um, which diagonalizes the effective mass 
matrix are then found algebraically via a method analogous to that of Ref. |53| . By making the 
substitution TJ°''^ , where C/"'^ = V^I^Um, the expression for the event rate at the detector 

is equivalent to Eq. |^ 

We again investigate the case of one additional non-standard interaction, namely a left-chiral 
scalar interaction. The standard and non-standard mixing matrices, U^^ and JJ^^^^ respectively, 
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are parameterised by three different rotation angles, 0i2, ^13 and 623 for U and </)i2, 013 and (^23 
for [/■S'i'S'i -pj^g values of the standard mixing angles are taken to be 



from Ref. |52j . The numerical evaluation of NMrec proceeds in the same manner as in Sec- 
tions [5I & IS] When considering muon production at the detector as a result of NSI we make the 
approximation that a^^^'^ ~ o"^^'^^. For the i^e ^ i^fj. channel the interference term in the cross 
section is ignored due to the comparatively small muon mass. 

In Figure^&Elwe have plotted the variation of NM^ec as a function of (f)i2 and ^13 respectively. 
Note that since we are only considering the effects of the electron background and, for the time 
being, no new physics at the detector the angle <l)23 decouples from the calculation as per the 
discussion of the previous section. As was expected the variation of NMrec in Figure 01 takes much 
the same form as in Figure |21 when — 0, or R{cf)i2) is diagonal. However when we choose 
4>i3 = 7r/4 and vary (j)i2 the results are dramatic. While not unexpected these results serve to 
highlight the point that any experimental analysis aiming to place bounds on the couplings of 
NSI must consider the full three neutrino system. This applies even for cases where the standard 
treatment may be done in a two neutrino framework. 

Also shown in Figure 01 are the effects of including non-standard physics at the detector in 
addition to the propagation stage. We see that the assertion made that the results of new physics 
at the detector are negligible when compared with effects of new physics during propagation are 
in this case justified. 

8 Concluding remarks 

In this paper we have united parameterisations of physics beyond the SM such as those of the 
Michel parameters and non-standard /3-decay couplings of precision weak physics experiments 
with the phenomenon of neutrino oscillations. As such a framework for analysing specific theories, 
such as Supersymmetry, leptoquarks etc. has been established. This was achieved by developing a 
scattering theory of the production, propagation and detection processes. In doing so we challenged 
the notion of what is termed a neutrino flavour state. This was necessitated by the fact that neutrino 
flavour states do not have a definite mass as is required of intermediate states in a scattering theory, 
and phenomcnologically there is no compelling reason to assume that new interactions couple to 
the same linear combination of mass states as the V-A interaction. As a result of this formalism we 
were able to show explicitly the known results that flavour violation or partial flavour violation can 
not lead to oscillation phenomena in a vacuum. And, in contradistinction, that flavour violation or 
partial flavour violation in matter can lead to oscillations even for degenerate vacuum masses. The 
condition for the latter result is that the non-standard interaction be non-diagonal in the standard 
basis. In addition we have shown that in the ultra-relativistic limit, production and detection 
processes with opposite chirality decouple, that is right-handed neutrinos become fully sterile from 
their left-handed counterparts. However in matter a right-handed potential acts as an effective 
mass linking the two states. 

The field theoretic description of neutrino oscillations was used to perform some simple calcu- 
lations for a generic neutrino factory. This type of experiment was chosen as a convenient starting 
point due the high intensity and energy of the beam with which the experiments are to be con- 
ducted. The beam energy is an important consideration since the field theory was derived in the 
limit that Sm^\L\/2Eq < 1. 

For the vacuum oscillation case we examined the — > Vr channel with a left chiral scalar 
coupling of strength G^^^^ = Q.OIG^^ . By allowing the non-standard mixing angle, 0, to vary we 
were able to observe flavour conserving and all flavour violating scenarios. We found the greatest 
effect, at high energies, was for flavour violating interactions which effectively picked out the non- 
standard term in the cross section. Just above the r-production threshold flavour conserving 
interactions play a more important role. 



sin 012 — 0.53 , 
sin 6*13 = 0.03 , 
sin 6*23 = 0.71 , 



(63) 
(64) 
(65) 
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The matter enhanced channel was also investigated. Two important differences from 

the vacuum case were observed. Although the coupling strength of the NSI was the same as that for 
the one used in the vacuum oscillation study the experimental sensitivity was up to three orders of 
magnitude greater, indicating that it is safe to ignore new physics at the detector and source for the 
matter enhanced channels. We also note that the greatest effect when varying the non-standard 
mixing angle, ij), did not correspond to direct flavour violation, rather when the non-standard 
interaction was maximally mixed in the standard basis. In the conventional parlance this would 
correspond to a decay to a superposition of flavour states. 

The focus of this paper has been one of developing a solid framework for future work. Future 
calculations will concentrate on more realistic experimental simulations than the ones presented 
here and on performing an analysis of the impact of a right-handed interaction on neutrino prop- 
agation through matter. 
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Figure 1: Top panel: The required number of detector mass-muon number units NM^-ec with 
qSlSl = O.OIG^^, dm? = 0.0025ey2, L = 732km. Bottom panel: The required number of 
detector mass-muon number units NMrec with G^^-^l = o.OlG^^, Sm^ = 0.00256^^, L = 732km. 
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Figure 2: Top panel: Relative contribution to the charged current event rate for the pure 
scalar term, long dash, the interference term, solid, and the sum of the two, short dash. With 
qSlSl = o.OlG^^, 6m^ = 0.0025eyVc4, L = Ti2km and = 50 GeV. Bottom panel: Relative 
contribution to the charged current event rate for the pure scalar term, long dash, the interference 
term, solid, and the sum of the two, short dash. With G'Si.Si, ^ o.OlG^^, 5m'^ = 0.0025eV^ /c^, 
L = 732fcm and = 20 GeV. 
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Figure 3: The required deteetor mass-muon number in units of 2 x lO^""^ kt/yr. The soHd Hne is 
with Ef^ — 50 GeV, the long-dashed line is with i?^ — 20 GeV. The vacuum mass difference is 
taken to be positive with \6m'^\ = 3.5 x lO^^eV^/c'', \L\ = 7332km and sin2(26') = 0.004. 




Figure 4: The required detector mass-muon number in units of 2 x 10^^ kt/yr and ioiE^ = 50 GeV, 
L = 7332 km, 0i2 = 0, 023 = ^23 and g^^^^ = O-Olg^^. The solid line is the event rate with no 
new physics at the detector. While the dashed line is the event rate with a non-standard coupling 
of G^^-s^ = O.OIG^^ at the detector. 
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Figure 5: The required detector mass-muon number in units of 2 x 10^^ kt/yr as a function 
mixing angle (pi2- The other parameters are -E^ = 50 GeV, L = 7332 km, 0i3 = 7r/4, 4)23 = 
and gr'St'Si, _ Q Q\gLL ^y^s plot we have assumed no new physics at the detector. 
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